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Condensation occurs in nonequilibrium steady states when a finite fraction of 
particles in the system occupies a single lattice site. We study condensation transi- 
tions in a one-dimensional zero-range process with a single defect site. The system is 
analysed in the grand canonical and canonical ensembles and the two are contrasted. 
Two distinct condensation mechanisms are found in the grand canonical ensemble. 
Discrepancies between the infinite and large but finite systems' particle current ver- 
sus particle density diagrams are investigated and an explanation for how the finite 
current goes above a maximum value predicted for infinite systems is found in the 
canonical ensemble. 



I. INTRODUCTION 

Condensation transitions fl| have been observed in a wide variety of nonequilibrium steady 
states . Such a transition occurs when a single lattice site contains a finite fraction of 
particles of the system. Condensation transitions are also manifested in a different guise in 
driven diffusive systems, where typically the particles obey hard-core exclusion. There the 
condensation appears as a transition from a freely flowing fluid phase to a jammed phase 
P]. Such jamming transitions have been reported, for example, in the context of traffic flow 
modelling m. 

A simple lattice model known as the Zero-Range Process (ZRP) p] illustrates an equiva- 
lence between condensation and jamming transitions and provides a framework within which 



condensation transitions may be analysed . In the model particles hop between lat- 

tice sites with rates which depend on the number of particles at the departure site. If the 
rates decay to zero or if they decay slowly enough to a finite value as the particle num- 
ber is increased, one finds a condensation transition whereby a single lattice site becomes 
macroscopically occupied. 

The connection between condensation and jamming transitions was first made in the 
context of the Bus Route Model ^ where an approximate description assigns hopping rates 
to buses that decay as a function of the distance to the next bus ahead. If the decay is 
suitably slow a jamming transition arises whereby the buses cluster together into a jam and 
there is a large region of the bus route devoid of buses. This approximate description of the 
bus route model maps exactly onto a ZRP where the buses correspond to the sites of the 
ZRP and the number of bus stops in front of a bus corresponds to the number of particles 
at a site of the ZRP. Phase separation in driven diffusive systems may also be related to 
condensation transitions and the ZRP has been used to formulate a general criterion for 
phase separation in one dimension 

In addition to the condensation mechanism in the ZRP described above (condensation 
associated with hopping rates that are decreasing functions of the occupation of a site) one 
can have condensation associated with disorder N, M . Here disorder refers to different sites 



having different hopping rates. The latter mechanism has been shown to be equivalent to 
Bose- Einstein condensation A particularly simple system that exhibits Bose-Einstein- 

like condensation is where there is one slow defect site at which a condensate forms at high 
enough density. 

In the present work we study a simple model, with a single defect site, that exhibits 
features of both condensation mechanisms (bus-route-like and disorder-induced) mentioned 
above. The steady state can be solved and the critical density for condensation in the 
thermodynamic limit can be predicted. However, finite systems show significant deviations 
from these predictions. In particular, the low density, fluid phase appears to continue to 
higher densities than would be allowed in the thermodynamic limit. This results in an 
overshoot in the current, above the thermodynamic saturation value. Related finite size 
effects have been reported in traffic flow models 

Our aim in this paper is to analyse in detail how the finite size overshoot in the current 



occurs. Now, since the ZRP dynamics conserves particle number the natural ensemble to 
consider is the canonical ensemble. However, one finds that analysis of the model is more 
direct within the grand canonical ensemble (i. e. where the particle number fiuctuates). It has 
been shown that in the infinite system limit the two ensembles are equivalent 0. However, 
since we are interested in the behaviour of large but finite systems, we broach the issue of 
the inequivalence of the two ensembles on finite systems. We are able to analyse the single 
defect site model in both the canonical and grand canonical ensembles and find that they 
give different predictions for the finite size overshoot in the current. In particular, within 
the canonical ensemble the overshoot is predicted to peak at a density y9 = + 0(M~^/^), 
where is the thermodynamic critical density and M is the number of sites, and we confirm 
this by simulations. 

The paper is organised as follows: in section 2 we define the ZRP and discuss the two 
types of condensation mechanism; in section 3 we introduce the single defect site (SDS) model 
and analyse the condensation transition within the grand canonical ensemble; in section 4 
we present exact numerical results and simulations that illustrate the overshoot of the low 
density fiuid phase on a finite system in the canonical ensemble; in section 5 we analyse the 
SDS model in the canonical ensemble; and we conclude in section 6. 

II. THE ZERO-RANGE PROCESS 
A. Model definition 

We consider a one-dimensional, asymmetric zero-range process. We have a one- 
dimensional lattice of M sites labelled ^ = 1 . . . M, upon which reside L indistinguishable 
particles, each site being able to hold any number (0 . . . L) of particles. Particles move from 
site to site under a dynamics which conserves the number of particles and is described as 
follows: a particle from site /z hops to site n + 1 with a rate, u^(n^), which can depend 
on both the departure site, /x, and the number of particles on this site, n^. The boundary 
conditions are periodic; a particle hopping from site M will hop to site 1. The configuration 
of the system at any given time is specified by the set of occupancies of each site, {n^}. 

Thus we consider the zero-range process on a finite system of size M; in particular we 



shall be interested in the limit of large M where the particle density if = L/M is held fixed. 



B. Steady state 

The one-dimensional zero-range process, with particles hopping to adjacent sites in one 
direction only, has a steady state which is straightforwardly solved 0| • Requiring that 
there be no net fiow of probability to or from a configuration, one finds the probability of a 
configuration, {n^}, in the steady state to be 

1 ^ 

where the f^ are given by 

I n -n-^ if n > 1 

W = ""'^""^ " (2) 

I 1 if n = , 

and Z(M, L) has been introduced to ensure the probabilities are correctly normalised. This 
normalisation is analogous to the canonical partition function from equilibrium statistical 
mechanics, being the sum of the un-normalised probability over all allowed configurations, 

(M \ M 

n^.K)- (3) 
fi=i / ^=1 

Here we have summed over all occupancies for each site, from zero to infinity, and included 
the delta function to ensure that only those configurations with the correct total number of 
particles contribute. Thus Z{M, L) is the normalisation in the canonical ensemble. 

This steady state probability distribution can be used to calculate other steady state 
properties of the system. For our purposes the most important of these is the average 
hopping rate from a site, v, which we refer to as the current. In the steady state this is 
independent of site and is given by 

- ^ E "(«o E ^ (t % - '] n - . w 

ni n2,...,nM \M=1 / M=l 

where in the final step we have used the fact that u^(0) = (the departure rate from an 
empty site is zero) and u^(x)f^(x) = f^{x — 1) for x > 1. 



C. Condensation transition - grand canonical analysis 



Depending on the chosen hopping rates the system can exhibit a condensation transition 
in the thermodynamic hmit M — > oo with the density of particles fixed. In this case, at 
low density the system is in a fluid phase, by which we mean that all sites have average 
occupancies that are infinitesimal fractions of the total number of particles. On the other 
hand, at high density the system is in a condensed phase wherein a single site holds a 
finite fraction of the total number of particles. We shall show that there are two distinct 
mechanisms by which this may occur. 

The simplest approach to analysing the condensation transition is to use the grand canon- 
ical ensemble. The functions ffiin) given in are only defined up to some multiplicative 
factor, z^, which we have tacitly set to one. Reinstating z, we can interpret z as the fugacity. 
Then we can get the equivalent of the grand canonical partition function in the usual way 



12j |. by summing over all configurations, including those that violate the constraint on the 



total particle number, and then choose z such that the average total particle number is equal 
to L: 

oo oo oo M 

ziM,z) = E E ■■■ E U'^'^M (5) 

ni=0n2=0 7iM=0 (1=1 
M 



U^.i^) (6) 



fi=i 

oo 

where F^(^) = J^z^f^in) . (7) 

n=0 

The constraint on the average particle number is written as follows 

M 

ij,=i 

with the average occupation of each site, (n^), being found as usual by: 

(«,) = . (9) 

which yields 

z J^F'Jz) 



ii=i 

where 9? is the particle density L/M. 



1. Different types of condensation 



To see the two apparently different mechanisms we analyse the behaviour of the rhs of 
(finjl as z ^ 2;max where Zmax is the radius of convergence of ((Zj). As z — > ^max one of two 
things can happen, either the rhs of (fTn|l will converge, or it will diverge. 

If the rhs of (fTUj) converges at z = -Zmax, then this will give a critical density, v^c, above 
which fllO|) can no longer be satisfied. Hence there is a phase transition, whereby the excess 
density condenses onto a single site. We refer to this as mechanism A. 

If the rhs of fllO|) diverges at 2; = ^max, then it is still possible to have a phase transition 
by a different mechanism. Recall that each term in the sum in (fTIH) is the density of particles 
at a site. It is possible for the rhs of (fTIHl to diverge due to only a single term, say site 1, 
of the sum diverging. This would imply that site 1 contains a finite fraction of the total 
number of particles, i. e. we have a condensate. We refer to this as mechanism B. To see 
that mechanism B is reminiscent of the Bose-Einstein condensation mechanism, one thinks 
of the particles as bosons and the sites of the ZRP as Bose states. Then site 1 into which 
condensation may occur corresponds to the ground state of the bosonic system. 

The choice of hopping rates will determine which mechanism takes place. Previously 
studied cases are: 

1. The hopping rates are the same for each site and decay with increasing particle number 

; then we can have a condensation through mechanism A. 

2. The hopping rates do not depend on the number of particles at a site, but do depend 
on the site, i. e. the system is disordered Then we can have a condensation through 
mechanism B. One very simple realisation of this is a single defect site with a hopping 
rate less than the other sites. Then at high density this site will support a condensate. 

One important difference between mechanisms A and B is that in A there can be a 
spontaneous symmetry breaking, in that the site supporting the condensate is selected at 
random. With mechanism B, on the other hand, the term in the sum which diverges specifies 
which site supports the condensate (taken as site 1 above). 

If we have a mixture of decay and disorder in the hopping rates, then it is possible to 
have a phase transition through either mechanism. In the present work we investigate this 



possibility. 



2. Overshoot in the current v 



In the grand canonical ensemble the current is given by 

V = {ui) 

^ oo oo M 

' ni=0 n2,n3,...,nj\/=0 /i=2 



where to proceed to the final line we have used the fact that u^(0) = and u^(x)f^(x) = 
i^{x — 1) for X > 1. 

Thus in the grand canonical ensemble we have an equivalence between the current, f, 
and the fugacity, z. From the discussion of the behaviour of z in the previous subsection we 
deduce that, in the infinite system limit, v increases with density in the fiuid phase until at 
density ipc the current saturates at f = ^max and we enter the condensed phase. 

However, on a large but finite system, things are more subtle. In mechanism B one can 
always solve ()10j) for z < ^^max- Thus v should always be less than its saturation value 
Zmax- In mechanism A, on the other hand, one cannot satis fy (TUj) for ip > ip^ with z < 
-2max- Thus it is possiblc that v overshoots the value 2;niax [13 1. In simulations (in the 
canonical ensemble) of finite systems an overshoot in the current versus density diagram 
can be observed, i. e. the current appears to go above the expected maximum value before 
dropping back down at higher density (see section HVjl . Thus on a finite system there can 
appear to be a continuation of the fiuid phase into the region where, on the infinite system, 
one would expect a condensate. 



III. SINGLE DEFECT SITE 



We now restrict ourselves to a very simple ZRP with all hopping rates constant and equal, 
except for a single slow defect site which has a hopping rate that decreases for increasing 
particle number. We refer to this as the single defect site (SDS) model. We begin by 
analysing the model within the grand canonical ensemble. 



We consider the one-dimensional lattice of /i = 1 . . . M sites with periodic boundary 
conditions from section III Al and now choose the following hopping rates (for I > 0) 

m{l)=p(i + ^] (12) 



Uf^il) = 1 for /i > 1 , (13) 

with p < 1 and we generally take p < 1/(1 + A), giving site 1 a hopping rate slower than all 
the others for occupancies / > 1. Then from 

1 ra' 

Un) = — — (14) 

f^(n) = 1 for /i > 1 , (15) 
where (a)„ is the Pochhammer symbol defined by 

{a)m = fl(o. + l)(a + 2) . . . (a + m — 1) 

(a)o = 1 . (16) 

Also, from ((7j) we have 

FiW = E(^) 7Y^ = 2Fi(l,l;l + A;^/p) (17) 

n=0 ^-^^ ^ 

F^(z) = F(;2) = for /i > 1 , (18) 

where the hypergeometric function is defined as 

OO / X 



2Fi(a,6;c;x) = > — r . (19) 

•i— ^ r ml 



rra=0 

To see the equivalence in (fTTjl note that (l)m = w^'- 

If we look at the grand canonical analysis, specifically equation (ITUI). for the SDS model 
the slow site gives a ratio of hypergeometric functions, while the 'non-slow' sites, /i = 2 . . . M, 
give a ratio of two simple geometric sums 



M VFi(^) ' ' F{z) 



Mp(l -f A) 2Fi(l, 1 ; 1 + A ; z/p) \ Mj 1-z 
<"'<^'> + fl-i,V„W>. (21) 



M V M 



In the last line {ni{z)) means the average number of particles at site 1, whereas {n{z)) means 
the average number of particles at a site /i 7^ 1. Thus, as M ^ 00, the critical density will 
be 

¥>c = ^ . (22) 
1 — p 

To see this we take the limit z z^i^^ = p and M ^ 00 of the expression ^IT\ such that 
(jii)/M 0. This last condition ensures we are not in the condensed phase. However, for 
Lf > ifcV^e must have {ni{z))/M finite to solve (j2^, thus we have a condensate on the slow 
site. 

In order to determine the finite size corrections to the current v = z we analyse (1211] for 

n 

M large but finite. We invoke the following results [1^ 

^ / , N T(c)T(c — a — b) , . 

if Oa + b hm 2F1 a,&;c;x = ^ \ ^. -( 23 

x^i- L [c — a) L [c — b) 

^ = " + ^ J-T- |ln(l-a;)| = mW) ^ ^ 

if c<a + b \nn p^l^ = ^Ml(-±]^ (25) 

where T{x) is the usual Gamma function. Thus 

if A>2 \na_{n,iz)) = (26) 

if A = 2 hm-il^lM) =1 (27) 



2— >p 



ln(l — z/p)\ 



if 1<A<2 \im{n,{z)){l-z/p)'-'= .^^ ^^'^ (28) 
2-*p- sm(^7r(^A — i J J 

if A = 1 lim_(ni(z))(l - z/p)\ ln(l - z/p)| = 1 (29) 

if A < 1 lim {ni{z)){l - z/p) = (1 - A) , (30) 



where to obtain ()28j) we have used the identity 



Therefore when — p) for A < 2 one can satisfy (PT|) by choosing 



if A = 2 I- z/p^ exp [-M((^ - v^c)] (32) 



if 1 < A < 2 I- z/p 



(A - l)2vr 



_sin(7r(A- l))(^-yPc)M 



l/{2-A) 



(33) 



if A = l ^ TT^^ (34) 

if A<1 l-z/p^^^—^^^. (35) 

For A > 2 on the other hand, we cannot satisfy ^21\ a.s z ^ p~ when > ^pc- Instead 
we must consider z > p. Thus the grand canonical treatment imphes an overshoot in the 
current v (= z) for A > 2. 

In order to consider z > p we impose cut-offs = L m the sums of (0)- This explicitly 
ensures that no site contains more particles than the total number L. Then 

(n.jz))^ - .fS-)(2,2;2 + A;./p) 
p(l + A) ,Yf\l^l-l + X-z/p) 

where ^F^ ^(a, 6 ; c ; x) = V ^^^^ . (37) 



m=0 ^ 



For A > 2 we write z/p 'm the following form: 

- = exp(a(L)). (38) 
p 

To satisfy (j2H) for ip > (pcwe require {ni{z)) = 0{L). It turns out that as z increases above 
p this is first achieved when the numerator of is 0{L) and the denominator 0{1). In 
anticipation of this we write 

(n,) ^ where G, = !!^E(^ dfc (39) 

^ ' A + Gx Ji 

and A = 0{1), i. e. we have replaced the sums in ()3(')P37|1 by integrals over the asymptotic 
forms of the summands, retaining a finite error term A in the denominator. Then we require 
that Ga = 0{1) and Ga-i = 0{L). Now the asymptotic behaviour of G^ is 

G,~^^g^ for La(L)»l. (40) 



Thus we require 



exp(La;(L)) 



0(l) 



(41) 



which imphes 



Lj 



InL Infln 




+ 0{l/L). 



(42) 



We can now summarise the resuhs of the grand canonical analysis for a large but finite 
system 

A For A > 2, (jii) remains finite as z — > Zmax and we identify this case with condensation 
mechanism A of section III C 11 In the condensed phase, in order to satisfy ()2H) z = v 
overshoots its saturation value p. The approach of v to its saturation value as particle 
number L increases is f/p ~ 1 + (A — 1) In L/L. 

B For A < 2, we can access the condensed phase through (j2T| and v approaches its 
saturation value from below as detailed in (jHSHHSl)- We identify the case A < 2 with 
mechanism B. 

IV. SIMULATIONS AND EXACT NUMERICS IN THE CANONICAL 

ENSEMBLE 

In this section we will present results obtained from Monte Carlo simulations of the 
system and exact numerical calculations of Z(M, L) in the canonical ensemble (fixed particle 
number). What we shall find is that on finite systems the overshoot in the current, v, is 
clearly present whenever A > 0. This is in contradiction to the grand canonical analysis of 
the previous section which predicts an overshoot only for A > 2. 

To observe the overshoot in the current, v, exact numerical calculations of Z(M, L) for 
systems of sizes up to M = 1000, with up to L = 2000 particles were undertaken. These 
calculations are performed recursively :2]. It is then possible to obtain v via the relation 

The current, v, is plotted against particle density at fixed system size for various values 
of A and p = 0.2 in Figured In this figure we see that for low density the current increases 
whereas for high density it saturates. For A > 0, in between these two regimes is an overshoot 
in the current where we observe non-monotonic behaviour. Remember that, for an infinite 
system, as the density is increased we would expect to see v increase to a maximum value 



and stay there, yielding a current versus density curve that is non-analytic at the critical 
density. Thus, for the infinite system, when v is increasing the system is in a fiuid phase, and 
when V saturates the system is in a condensed phase. From Figure Q we see that on a finite 
system, for A > 0, f overshoots above its expected saturation value. This effect increases 
with increasing A. 

Notice also that for a finite system, v crosses its predicted saturation value for an infinite 
system, in this case p, close to the predicted critical density for an infinite system, given by 
from equation (j221)- The current, f , is also shown in FigureEl plotted against particle 
density, but for different system sizes, M = 200, 400, . . . , 1000 with A = 4 and p = 0.2, from 
which it can clearly be seen that the overshoot decreases in severity for increasing system 
size. Thus the current versus density curve will approach the expected infinite system result. 

0.25 I ' , ' , ' , ' 1 
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FIG. 1: Exact numerical calculation of the current, v, plotted against particle density for a system 
of size M = 1000, with p = 0.2 and A = 0, 1, 2, 3, 4 from bottom to top. 



In order to investigate what is actually happening in the overshoot region, Monte Carlo 
simulations of the system were run at chosen points within the overshoot. The simulations 
were run in three key areas: firstly, between the predicted critical density on an infinite 
system (j22|) and the maximum of the overshoot; secondly, between the maximum of the 
overshoot and where v drops down close to its expected maximum value of p; and thirdly 
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FIG. 2: Exact numerical calculation of the current, v, plotted against particle density for a system 
with A = 4, p = 0.2 and varying system size M = 200, 400, . . . , 1000 from top to bottom. 

where v has dropped down close to p. As we are interested in condensation- like phenomena, 
the simulations output graphs of the number of particles on the slow site with time. Some 
typical examples of these graphs are shown in Figure El 

Of course, on a finite system, strictly one cannot identify thermodynamic phases. How- 
ever, it is useful to think of putative phases emerging in large systems. In the following and 
Figure 01 'phase' is often taken to mean such a putative phase. 

Inspecting Figure IHl we see that the overshoot is initially an extension of the fluid phase 
(Figure Et); up to the maximum of v the system is still in the fluid phase. At the maximum 
a condensate begins to emerge on the slow site. In the range of densities just beyond the 
maximum of the overshoot the emergent condensate appears intermittently, but is not stable 
over long time periods (Figure Eb)- We interpret this as 'coexistence' between the fluid and 
condensed phases, i. e. at any given time the system is in one or the other. Eventually this 
coexistence gives way to a clear condensate as v drops back down to its predicted maximum 
value of p (Figure Et)- 

Having established what takes place in the overshoot region, we proceed in the next 
section to give an analytical account for this effect working within the canonical ensemble. 
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FIG. 3: Graphs showing the number of particles on the slow site against time for A = 4, p = 0.2 
with M = 1000, obtained via Monte Carlo simulation, (a) System is at (/9 = 0.29, between expected 
critical density, in this case 0.25, and maximum of the overshoot. This is an extension of the fluid 
phase, (b) System is at = 0.34, just past the maximum of the overshoot. This is the coexistence 
region, a small condensate is seen to appear intermittently, but is not stable over long periods of 
time, (c) System is at (/? = 0.40, beyond the 'end' of the overshoot region where the current is 
approaching its expected maximum value v = p. Here the large pile of particles at the slow site is 
stable and the system is in the condensed phase. 



V. CANONICAL ANALYSIS 

The SDS detailed in (fT^ [T^ . has a convenient property in that the canonical normali- 
sation, Z(M, L) 0, can be reduced to a single sum. This makes it possible to analyse the 
model directly within the canonical ensemble. 

The form that Z(M, L) takes for the SDS model is 

ZiM,L) = ±p-^(^^ "-"-'). (43) 



Here we have simply inserted the ()14I15|) into the general expression for Z(M, L) Q. The 
simple form of the hopping rates for sites fi > 1 (fT^ allows the sum over n2. . .um to be 
performed easily, yielding the combinatoric factor which comes from the delta function in 
(jH))- We also re-label ni as n for convenience. 

The normalisation is the sum over all un-normalised probabilities of occupancies of the 
sites. This means that the magnitude of the n^^ term in the sum ()43|) is proportional to 
the probability of the slow site having occupancy n. Thus if Z{M,L) is dominated by a 
maximum at low n, the corresponding system is in the fluid phase. Conversely, if Z(M, L) is 
dominated by a maximum at large n, the corresponding system is in the condensed phase. 
If there are two maxima of similar magnitudes we interpret this as phase coexistence. Thus 
investigating the maxima of this sum should furnish us with an explanation of the overshoot. 

To investigate the maxima of this sum, we look for its stationary points, i. e. points where 
the ratio of two consecutive terms tends to one. Solving for these points gives a quadratic 
equation in n, the solutions of which are: 



M , ^ A M 



, (44) 

where = p/ (1 —p) is the expected critical density and the constant A = ipc/p + (fd^ + A). 

The nature of these turning points can be determined by considering the ratio of the first 
two terms and also of the last two in ()43|) . In the following we also assume that p < 1/(1 + A) 
for simplicity, although this turns out to make little difference to the overall picture. 
First we note that the roots in (jl^ become real at a density ip2 given by 

AXp 1 8Xp 



p 1 /l+p(l-A) 

<^2 - 



1-p M \ 1-p 



M{i-py M^{i-py 



(45) 



Thus on increasing the density from zero the profile of the terms in the sum goes through 
the sequence illustrated in Figure IH (a) For < ip2 there is a boundary maximum at n = 0; 
(b) at = a stationary point emerges; (c) for (f > maximum emerges at large 

n in addition to the local maximum at = 0; (d) the maximum at large n increases and 
dominates, with n = eventually becoming a boundary minimum. 

We interpret this sequence as follows. At low density, (f < (f2 there are no turning points 
in the physical region < n < L, with n real. Thus the values of n near n = dominate and 
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FIG. 4: Graphs showing the logarithm of the terms from the normahsation sum Z(M, L) ()43|). for 
a system of size M = 1000, with A = 4 and p = 0.15. (a) L = 200 {ip < ip2), boundary maximum 
at n = 0, no turning points, (b) L = 233 {ip = (p2), boundary maximum at n = and stationary 
point at higher n. (c) L = 300 {ip > ip2), boundary maximum at n = with greater maximum at 
high n. (d) L = 500 (99 ^ ip2), boundary minimum at n = and maximum at high n. 

this corresponds to the fluid phase. If we were to consider the M —>■ 00 limit of ()45p we would 
obtain (f2 = fc- However for M finite, the square root part of ()44|) remains imaginary up to 
(p2 > v^c and the system remains in the fluid phase until this point. Thus the maximum at 
high n emerges beyond the expected critical density = p/ (1 —p) by an amount (9(M~^/^): 

^2 = ^ + ^f^V'%0(-). (46) 

1-p i-p\mJ ^m' ^ ' 

When ip is just greater than Lp2 we interpret the boundary maximum at n = and the 
turning point at large n as coexisting fluid and condensed phases. Finally, for Lp much larger 
than (p2 the maximum at high n dominates and we interpret this as the condensed phase. 



We expect that the maximum of the overshoot will be near to the point as this is where 
a putative condensate can begin to form and reduce the current. We see that an amount 
(9(M~^/^) beyond is consistent with the data from exact numerical calculation of the 
current (Figure E]). In fact, we compared the peak of the overshoot with ip2 for systems of 
size M = 1000 we found very good agreement for A ~ 1 ... 6. Outwith this range we believe 
that further finite size effects are becoming important. 

The form of Z(M, L), given in equation (pSj) . also allows the behaviour of the current in 
the condensed phase to be analysed. One finds that the asymptotic behaviour of the current 
for large particle number (i. e. where the sum in is dominated by a maximum at large 
n) is 



P 



1 A /I 



(47) 



Thus for A > the current approaches its asymptotic value from above as oo. It should 
be noted that the second term in the expansion is small only for — ip^^ X/M. 

VI. DISCUSSION 

In this work we have studied the condensation transitions of a ZRP with a single defect 
site, in both the grand canonical and canonical ensembles. Analysis in the grand canonical 
ensemble predicts the two different condensation mechanisms (A and B) discussed in Sec- 
tion III C 11 For A > 2 mechanism A applies which predicts a finite size overshoot in the 
current. The approach of v to its saturation value p is v/p = 1 + 0{lnL/L). For A < 2 
mechanism B applies and there should be no overshoot. 

However, simulations in the canonical ensemble reveal an overshoot for all A > 0; only 
for A = do we see the expected mechanism B behaviour. Analysis within the canonical 
ensemble confirms this and predicts that for A > the overshoot in the current peaks at a 
density = + 0(M^^/^). The approach of v to its saturation value as particle number 
increases is v/p = 1 + 0{1/L). Thus it is not clear to what extent the mechanisms A and 
B are relevant or well-defined within the canonical ensemble. Moreover, from the point of 
view of finite size critical behaviour, the two ensembles are not equivalent. 

It would be interesting to further understand the finite size inequivalence of the ensembles. 



In order to shed some light on this one can follow a standard approach 







I2L Il5l | to evaluating 



the canonical normalisation Q in terms of the grand canonical partition function (0). We 
use the integral representation of the delta function to write the canonical normalisation as 
a contour integral in the complex plane 



Then for large system size and particle number, M and L respectively, this integral will be 
dominated by its saddle point, given by 



where is the particle density (L/M). This equation precisely recovers ()10|) . However it 
must be borne in mind that P9|l holds only when the saddle point of P8|) is valid, i. e. when 
it dominates the integral ()48|) . It would not be surprising for the saddle point to break down 
when it approaches the radius of convergence of Z{s) in a way that depends on system size, 
as in the analysis of ()10|) in section III CI Further analysis of this phenomenon would be 
enlightening. 

Finally it would be of interest to perform simulations within the grand canonical ensemble 
to confirm our analytical predictions. To do this one would introduce transitions whereby 
particles are created and annihilated with appropriate rates. 
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